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® Radiation passing through a cone beam collimator is received by a radiation detector (10) such as a.gamma 
camera head, as the gamma camera head Is moved In a helical orbit. Data g(n.u,v) collected during the helical 
orbit is scaled (42) to scaled data G(n,u,v). A first partial derivative dG(n,u,v)/6u is taken (46u) with respect to a 
horizontal direction and a second partial derivative dG(n,u,v)/6v Is taken {46v) with respect to a vertical direction. 
The partial derivatives are linearly combined (48) by being multiplied by sine and cosine values of an angle a 
between the horizontal direction u and an arbitrary direction p in the detector plane to form partial derivatives 
dG(n,u,v)/dp. The coordinate system of the derivatives is converted (60) from the (n,u,v) coordinate system to an 
(n,a,p) coordinate system. The first derivatives are projected (62), I.e. summed row by row, onto a q axis which 
is perpendicular to the p arbitrary direction. The one-dimensional projection arrays are rebinned (70) to form first 
derivative Radon domain data R*(0.p). A second derivative R"(9, p) Is taken (74) of the Radon domain data. The 
second derivative Radon domain data is backprojected (80) into an image memory (82) and displayed on a 
video monitor (84). 
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This invention relates to imaging methods and apparatus. More particularly the invention relates to an 
apparatus and method for generating an image representation of an Interior portion of a subject. It finds 
particular application in conjunction with single photon emission computed tomography (SPEC!) scanners 
with cone-beam collimation for medical diagnostic imaging and will be described with particular reference 

5 thereto. It is to be appreciated, however, that the invention will have other applications in which cone-beam 
type data is reconstructed into an image representation for medical, quality assurance, and other examina- 
tions. Although described in conjunction with emission radiation sources which emit radiation from the 
subject, it will be appreciated that the present invention is also applicable to reconstructing image 
representation from transmission radiation sources which transmit radiation through a subject. 

10 Cone-beam collimators are commonly used with single photon emission computed tomography and 
other gamma camera devices. The cone-beam collimator diverges outward from a subject toward the 
scintillation crystal or detector head. This enables a large fraction of the detector head crystal face to be 
utilized when Imaging relatively small regions of the patient, e.g. the heart. This effective magnification 
produces a combination of resolution and sensitivity gains over images formed using parallel or fan beam 

75 collimators. 

Most commonly, cone-beam projection data is converted into image representations using a technique 
developed by Feldkamp, Davis, and Kress which is described in "Practical Cone-Beam Algorithm", J. Opt. 
Soc. Am. Vol. I, pp. 612-619 (1984). The Feldkamp technique uses an algorithm which was derived using 
approximations of a fan beam formula. A fan beam commonly lies in a single plane, but diverge in that 

20 plane. When multiple fan beams are used concurrently, the planes are substantially parallel. In this manner, 
the radiation paths diverge along one axis and are parallel in the other axis of the plane of reception. 

Feldkamp, et al. use a convolution and backprojection method which assumes that the focal point orbit 
Is a circle. However, If the foca! point of the collimator follows a single, planar orbit, the obtained data is not 
sufficient for an exact three-dimensional reconstruction. The insufficiency in the amount of collected data 

25 causes distortions and artifacts in the resultant image. There are three kinds of artifacts in the images 
reconstructed by the Feldkamp, et al. algorithm: a) reduced activity on non-central slices; b) cross-talk 
between non-central adjacent slices; and. c) undershoots in the transverse direction. 

In order to generate a complete or sufficient set of data, every plane which passes through the imaging 
field of view must also cut through the orbit of the focal point at least once. See Tuy "An Inversion Formula 

30 for Cone-Beam Reconstruction", SI AM J. Appl. Math. Vol. 43, pp. 546-552 (1983). The single planar orbit of 
Feldkamp does not satisfy this condition. 

Once the projection date is complete, one reconstruction approach is to convert the cone-beam 
projections to Radon transforms and use the Radon Inversion formula to reconstruct the image. This 
technique involves rebinning or sorting of the cone-beam data into another format. See Grangeat "Analysis 

35 d'un Systeme D'Imagerie 3D par Reconstruction a Partir De X en Geometrie Conique" Ph.D. Thesis I'Ecole 
Natlonale Superieure Des Telecommunications (1987). 

Others have proposed mathematical improvements to the reconstruction algorithms. For example, the 
cone-beam data sets can be inverted if one assumes that for any line that contains a vertex point and a 
reconstruction point, there is an integer M (which remains constant for the line) such that almost every 

40 plane that contains this line intersects the geometry exactly M times. See Smith "Cone-Beam Tomography: 
Recent Advances and a Tutorial Review", Optical Engineering, Vol. 29 (5), pp. 524-534 (1990). However, 
this integer requirement condition is too restrictive for practical application. 

According to a first aspect of the present invention there Is provided an apparatus for generating an 
image representation of an interior portion of a subject, the apparatus comprising: a radiation detection 

45 means for receiving radiation travelling along a cone of rays which converge at a focal point and for 
generating electrical data indicative thereof; a means for moving the radiation detecting means such that the 
focal point moves in a helical orbit and the radiation detection means generates helical orbit data; a partial 
derivative means for generating a partial derivative of the helical orbit data with respect to an arbitrary 
direction in a plane of the detector means; an integrating means for projecting the derivative along a 

50 direction perpendicular to the arbitrary direction in the detector plane; a rebinning means for rebinning the 
projections into first derivative Radon domain data; a second derivative means for generating a second 
■ derivative of the Radon domain data; a backprojecting means for reconstructing the second derivative 
Radon domain data Into a three-dimensional image representation; and an image memory means for storing 
the three-dimensional image representation. 

55 According to a second aspect of the present invention there is provided a method for generating an 
Image representation of an Interior portion of a subject, the method comprising: generating cone beam 
electrical data from detected radiation travelling along a cone of rays which converge at a foca! point, the 
radiation being detected by a radiation detection means along a detection plane; moving the radiation 
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detecting means such that the focal point nnoves in a helical orbit: taking a partial derivative of the electrical 
data with respect to an arbitrary direction along the detection plane; projecting the partial derivative along a 
direction perpendicular to the arbitrary direction along the detection planet; rebinning the projections Into 
first derivative Radon domain data; taking a second derivative of the Radon domain data; backprojecting the 

5 second derivative Radon domain data into a three-dimensional image representation; converting at least a 
portion of the image repesentation into a human readable display. 

One advantage of an apparatus and method embodying the present invention is that the resultant 
images are reconstruction algorithm artifact free. 

Another advantage of an apparatus and method embodying the present invention is that the reconstruc- 

10 tlon is exact. 

Another advantage of an apparatus and method embodying the present invention is that the raw data 
can be generated using conventional hardware which moves detector heads in a circular orbit and linearly 
along the patient. 

An apparatus and method in accordance with the present invention will now be described, by way of 
75 example, with reference to the accompanying drawings in which : 

Figure 1 is a diagrammatic illustration of the apparatus which is SPECT imaging system; 
Figure 2 illustrates the geometry of a helical scan by the apparatus; 
Figure 3 illustrates a cone beam projection in use of the apparatus; and 
Figure 4 illustrates a Radon transform In use of the apparatus. 
20 With reference to Figure 1 . a cone-beam radiation detection means, such as one or more gamma 
camera heads 10, each with a cone-beam collimator 12. is mounted to move rotationally around and linearly 
along an examination region 14. The cone-beam collimator 12 has a plurality of paths defined by bores in a 
solid lead sheet or by lead vanes which focus at a focal point 16. The cone-beam collimator is oriented 
such that the focal point 16 is disposed across the examination region 14 from a subject therein. Radiation 
25 emanating from the subject or passing through the subject follows diverging paths through the collimator 12 
to the gamma camera head 10 or other detector. In this manner, a relatively small region of the subject is 
projected onto a relatively large region of a crystal face of the detector head 10, i.e. an effective 
magnification. 

The detector heads 10 are mounted on a gantry means or portion 20. The gantry means includes a 

30 plurality of motor drives 22 which can be operated individually or in combination in order to move the 
detector heads along selectable orbits. The heads are rotated along a circular orbit and the gantry is 
translated to move the heads axially along the subject An orbit controller 24 generates motor control 
signals for each of the motors to cause the beads to move along the selected orbit. More specific to the 
illustrated embodiment, the orbit controller includes a look-up table 26 which is preprogrammed with the 

35 appropriate positions which the detector heads and gantry should take to move the focal point 16 along the 
spiral orbit. A current position sensing means 28 monitors the current position of the detector head(s), such 
as by monitoring the angular position around the subject, the radial position toward and away from the 
subject, and the longitudinal position along the subject. A comparing means 30 compares the look-up table 
values with the actual rotational and longitudinal positions of the detector heads and gantry. A series of 

40 motor drivers 32 supply motive power to the motors or linear drive means 22 until the monitored current 
position matches the desired position from the look-up table. Optionally, descriptions of a plurality of orbits 
are selectively loadable into the look-up table 24. 

The gantry or an associated control console includes a data processing means for processing the 
output data from the detector head(s). More specifically, each detector head conventionally includes a 

45 scintillation crystal that is viewed by an array of photomultiplier tubes. Each time a radiation event occurs, 
the radiation passing through the collimator and striking the crystal causes a light flash or scintillation. The 
photomultiplier tubes nearest the scintillation respond with proportional output signals. Position and energy 
resolving circuitry connected to the photomultiplier tubes detenfnine the energy and position, hence the ray 
or direction along which the radiation travelled from the radiation event within the subject through the 

50 collimator to the detector head. Due to the cone-beam collimator, there is a direct relationship between the 
position on the scintillation crystal at which the radiation was received and the directional vector n of the 
ray. See FIGURE 2. 

Output data g(n,u,v) from the detector head during the helical orbit r is stored In a projection data 
memory 40. In the projection data coordinate system, n is the index of the views, i.e. an angular orientation 
55 around the helix; v is an axis along the patient, i.e. parallel to the patient support couch; and u is coordinate 
of the detector head 10 perpendicular to the axis v. In the preferred embodiment, the axes u and v are 
horizontal and vertical axes, respectively. A scaling means 42 scales or weights the projection data to 
generate weighted projection data G(n,u,v) which is stored in a weighted projection data memory means 
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44. Preferably, this scaling or weighting means multiplies the projection data g(n,u,v) by a constant based 
on the geometry of the gantry and the system as described in Equation (9) below. 

A first partial derivative means 46u takes a partial derivative of the weighted projection date G(n,u,v) 
with respect to the direction of axis u. A second partial derivative means 46v takes the partial derivative of 

5 the weighted projection data G(n,u,v) with respect to the direction v. A linear combining means 48 linearly 
combines the partial derivatives with respect to u and v. More specifically, the linear combining means 
includes a cosine multiplying means 50 and a sine multiplying means 52 for multiplying the partial 
derivatives by the cosine and sine, respectively, of an angle a between the axis u and an arbitrary vector 
direction p. An adding means 54 adds the products of the partial derivatives of the weighted projection data 

10 memory and the sine and cosine values to form a partial derivative of the projection data along the p 
direction, 6G(n,u,v)dp. See Equation (17) below. With reference to FIGURE 3. this rotates the partial 
derivative within the plane (u,v) by a degrees to fall along an arbitrary vector direction p. The partial 
derivatives of the projection data along the vector direction p are stored in a partial derivative along the 
vector direction p memory means 56. Thus, the means 46-56 function as a means 58 for taking a partial 

15 derivative with respect to the vector direction p. 

For simplicity, the variables of the partial derivative with respect to p have been redefined dG(n,a,p)/dp. 
As explained above, u and v are related to the unit vector direction p by sine and cosine relationships of 
the angle a (see FIGURE 3). However, in the (n,a,p) coordinate system, the summation along the q vector 
direction Is more straightforward. A coordinate transform means 60 transforms or rotates the (n,u,v) 

20 coordinate system of the partial derivative to a coordinate system (n,a,p). 

An integrating means 62 integrates or projects the data along a vector direction q which lies in the 
(u,v) plane perpendicular to the vector direction p. The integration is preferably performed by summing row 
by row or line by line to obtain a one-dimensional array along each vector direction q. See Equations (10) 
and (11) below. The plurality of one-dimensional arrays generated at each angle a are stored In a projection 

25 date memory means 64. 

With reference to FIGURE 4, the relationship between the integrals or projection data and the first 
derivative of the Radon transform R*(?,p) is given by Equation (12) below, where /9 is the angle between 
lines S'C* and O'C. A rebinning means 70 reblns or sorts the two-dimensional data arrays from the 
projection memory means 64 Into the first derivative R'(7,p) of the Radon transform for storage in a first 

30 Radon transform derivative memory means 72. 

A differentiating means 74 takes the second derivative in the Radon domain to generate second 
derivative Radon domain data H"(e,p) which is stored in a second Radon domain derivative memory 
means 76. A backprojecting means 80 backprojects the second Radon domain data arrays R"{e,f>) into a 
three-dimensional image representation f(x,y,z) which is stored in a three-dimensional image memory 

35 means 82. Data from the three-dimensional image memory means is selectively withdrawn for display on a 
video monitor 84. As is conventional in the art, planes of data may be selected and displayed. Alternately, 
surface renderings, projections, and other manipulations of the data may be made prior to display, as are 
known In the art. Data from the image memory means 82 may also be stored to disk or other archive 
memory means, 

40 Looking to the theory behind the reconstruction and with reference again to FIGURE 2, the cone beam 
projections are defined as line integrals: 



g(8',a) « j£{fl^tH)dt (1), 



where S' is the focal point on the orbit r and is the unit vector of the integral line. 

Using the cross point A of the Integral line and the detector plane as a parameter, the cone beam 

50 projections can be written as: 



g(a^x) - /£(8'*ti)dt (2). 

55 — 

In FIGURE 2, ( is the detector plane, O'-uvw is the detector coordinate system, 0-XYZ is the 
coordinate system for the object f(x,y,z), Rfan is the focal length, and r is the orbit of the focal point. 
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With reference to FIGURE 4. a 3D Radon transform R(?,p) of the object f(x,y,z) for a definite 8 , and 
p is defined by the plane integral: 



P) • / fft (2) A (2 • 8-p) d2 (3) , 



10 where 



6 « (•inOeoB^fSl&esin^^cosO) (4) 

75 



is the unit vector of the normal direction of integral plane a{e ,p). 

The first derivative of the Radon transform with respect to p is given by: 

20 



tL^a.p} ' ^u9,p) (5). 



25 



From alt the Radon transforms of an object, one can reconstruct the object exactly by using the inverse 
Radon transform. The inverse Radon transfonm can be expressed as: 

35 where: 

dc«» = sin0d0d$ (7). 

The first derivative of the Radon transfom^ can also be used to express the Radon inversion formula: 



40 



t (« - / -^E'ce, Pi,*«)d» (8) . 



45 

Thus, from the first derivatives of the Radon transforms, one can recover the object exactly. Grangeat's 
formula relates the cone beam projections and the first derivative of the Radon transform. 

With reference to FIGURES 3, consider a source position S' on the orbit r and a given unit vector t In 
2ir. Rrst, the weighted cone beam projection G(S',A) is defined as: 

50 



OtS'.H) - -^g(8',X) (9) 



55 



The integration plane a(e,p) is defined as the unique plane perpendicular to the vector t and crossing S'. 
The point C. the orthogonal projection of the origin O on this plane, is called its characteristic point. 0O*p. 
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Line AA'C'BD Is the cross line of integration plane a and detector plane ^ Notice that C* is the orthogonal 
projection of the origin O' on the line AA'C'BD. Line O'C is perpendicular to the line AA'C'BD. Connecting 
S' and C gives S'C From O', draw O'C" perpendicular to S'C. It Is very easy to prove that O'C" Is 
parallel to the OC and both of them are perpendicular to the plane o. One can also show that the points C, 

5 C", and A' are in the same line. There are now three coordinate systems: 0-XY2 is the coordinate system 
for the object f(x,y,z) and its Radon transform, O'-uvw is the coordinate system for the focal point and the 
detector plane, and O'-pqw is the rotation of O'-uvw along axis w by an angle a to set axis p perpendicular 
to the line AA'C'BD. hence the integration plane. 

The SG(S',?) is defined as the integration of the weighted cone beam projection G(S',A) over the 

70 integration line AA'C'BD: 

80(8^0) - f 0(8^A)d& (10). 

75 

The coordinates system O'-pqw can also be used to express the integral above as: 



20 



80(8^0) « Jo(8^0C^*q3)dq (11) 



where q is the unit vector of the axis q. 
25 Grangeat's formula gives the fundamental relation between the cone beam projection and the first 
derivative of the Radon transform R'(?,p): 



30 



(12), 



35 



40 



45 



where ^ is the angle between line S'C and O'C, the left hand side of the Equation is in the cone beam 
projection domain, and the right hand side is in the Radon domain. 

It is known that If the focal point follows a planar orbit, then the data obtained does not satisfy Tuy's 
sufficiency condition for exact three-dimensional reconstruction. Tuy's sufficiency condition says: If any 
plane that intersects the object f(x,y,z) intersects the orbit of the focal point and the detector plane at the 
same time at least once there is enough cone beam projection data to recover f(x,y,z). A helical orbit of 
measurement fulfills Tuy's condition. 

For example, a helical orbit with a pitch of 32 is selected with 64 views of projection data generated for 
each pitch. Thus, two pitches provide 128 views of projections which are uniformly distributed on the orbit. 
In the present example, the 0-XYZ is the coordinate system of the object f(x,y,z) which now is the Defrise 
phantom. In the 0-XYZ system the Equation of the helical orbit r can be written as: 



50 



Rfan * C08Y 
Rfan * 8lnT 

T— - 32 



(13), 



where y is the parameter of the Equation. It is the angle between the orthogonal projection of the focal point 
55 to the 0-XY plane and the x axes, as shown in FIGURE 2. 

In the coordinate system 0-XYZ the Equation of the integral plane o(e,p) is expressed as: 
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20 



30 



35 



40 



2 • 6 » p (14) . 

5 By using Equation (4), Equation (14) can be written as: 
xsindcos</> + ysindsin^ + zcosd - p « 0 (15). 

The focal point positions 7 (views) are obtained by solving Equations (13) and (15) for given point (e ,p). 

70 With reference to FIGURE 1, the helical cone beam projection data g(S',A) is stored in the three- 
dimensional array g(n,u,v) in the data memory 40, where n from 1 to 128 is the view of projection which 
means the location of point S", u. and v are horizontal and vertical components of point A in the detector 
plane. The detector plane In the present example is a 64x64 two-dimensional array. 
The projection data g(n,u,v) is scaled by means 42 to G(n,u,v) as per Equation (9). 

75 In order to use Grangeat's formula (12), the 6G{n,u,v)/&p of G(n,u,v) is obtained. It is known that: 



m m 

^80(8^9)- ^/o{8^oc^*q3)dq-/-^o(a^oc^*«9)d« (16) 



where G(S',OC'+qq)can be expressed as G(n,u,v). 

First, the partial derivative with respect to the p is calculated by means 58 and then the integral along 
25 the q direction is calculated by means 62. 

For the given tilt, angle o of the axis p, the partial derivative in the direction p is given by: 

-^0(n,u«v) = co8«-^g(n,u,v) 4- 8lns-^0(n,u,T) (17), 

which is calculated by means 46u, 46v, and 48. The line integral along AA'C'BD 



can be evaluated by using the algorithm described by P.M. Joseph in "An Improved Algorithm For 
Reprojecting Rays Through Pixel Images". IEEE Trans, on Med. Imaging MI-1, (3) 192-196. After the 
calculation above, the three-dimensional array !:dG(n^,p)/bp is calculated 62 and stored in memory means 
64, where n is the index of the views, a e [0*. 180'] is the tilt angle of the axis p, p is the coordinate in the 
45 p axis of the detector coordinate system 0-pqw. i:dG(n,a,p)/6p is the first derivative along p axis of the 
projection at (n^,p). 

The rebinning means 70 converts the three-dimensional array rdG(n,a,p)/bp to the three-dimensional 
array R'(?,p) in the Radon domain which is the first derivative of Radon transform of f{x,y,z) at (?,p) in the 
coordinate 0-XYZ. 

60 With reference to FIGURE 3. from given point (?,p) and Rfan one can solve Equations (13) and (15) for 
S'. or at least one point of S*. This can be converted to n (the views of projection). Line OO' is the Z 
component of the point S'. Now one uses (f ,p) and OO* to express p and a. From FIGURE 3. it can be 
seen that: 

55 O'C"=p-OO'cos0 (18), 

and because of 0'C"lS'C": 
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O'C" _ Rfan /-av 



where O'C" is given by Equation (18). Similarly, the expression of angle a is: 



10 



(20) 



75 

Equations (19) and (20) are the rebinning equations which are performed by the rebinning means 70. The 
rebinning means 70 further solves Equations (13), (15) and uses Grangeat's formula to get the first 
derivative of Radon transform R'(?,p) that is stored in memory 72. 

The first rebinning technique finds (n,a,p) for given (e,p), i.e. maps to the cone beam domain from the 
20 Radon domain. An alternative rebinning technique finds (?,p) for given (n,a,p), i.e. maps to the Radon 
domain from the cone beam domain. From FIGURE 3, It is apparent that: 



O'C • Rfan = O'C" • VRfan)* + (0'C')2 (21). 
25 That is: 



/(R£an)> ♦ (O'CO* 

30 



From FIGURE 3, It is also apparent that: 

35 



_JE_ , >^'0' * 00' X 1 ♦ OOi (23) 



40 

Solving for p gives: 



p . o'c'd ♦ - oWl ♦ ^^^) (24) . 



By using Equation (22). Equation (24) can be written as: 

50 

p « (ofc' ♦ OO'flliE) (25) . 

/(Rfan)' ♦ (O'O' 

55 

Because O'C is just p and 00' is easy to get by a given n, one can get p by solving Equation (25). 
First, the angle $ is calculated: 



9 



EP 0 592 093 A2 



cose « ^ - • (26) . 



Solving for 9: 



70 

6 = aeosl 



/ R£an >sln» N ^27) 



75 The next step is to calculate the angle <>. Because the point S' is on the plane o of FIGURE 3: 

• 8 = p (28) . 

20 

Combining this Equation with Equation (4) and Equation (13): 

R£an • fllnS • cos(y-4) = P ♦ (32 -7—) • cos9 (29) . 

25 * 



Giving the expression for ^: 



30 



35 



(32-yM) 



- CO0O 



Rfan - slnO 



(30), 



where 7 is obtained from a given n. Equations (25), (26), and (30) form the second rebinning equation set, 
which enables one to obtain the first derivative of Radon transform R'Ce'ip) by using Grangeat's formula. 

From R'(0',p). the second derivative means 76 takes the second derivative and the backprojection 
means 78 backprojects the second derivative to get the object f(x,y,z). In summary, the reconstruction 

40 scheme is: 



45 



Because the helical orbit produces a complete set of data, an artifact-free image can be reconstructed. 
50 Claims 



1. An apparatus for generating an image representation of an interior portion of a subject, the apparatus 
comprising: 

a radiation detection means (10) for receiving radiation travelling along a cone of rays which 
55 converge at a focal point (16) and for generating electrical data indicative thereof; 

a means (20) for moving the radiation detecting means (10) such that the focal point (16) moves in 
a helical orbit and the radiation detection means (10) generates helical orbit data (G(n,u,v)); 

a partial derivative means (58, 60) for generating a partial derivative (dG (n. a, p)/dp) of the helical 
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orbit data (G(n,u,v)) with respect to an arbitrary direction (p) In a plane of the detector mans (10); 

an integrating means (62, 64) for projecting the derivative {6G (n, a, p)/dp) along a direction (q) 
perpendicular to the arbitrary direction (p) In the detector plane; 

a rebinning means (70. 72) for rebinning the projections into first derivative Radon domain data (R' 
5 (e.p)); 

a second derivative means (74, 76) for generating a second derivative (R"(?,p)) of the Radon 
domain data; _^ 

a backprojecting means (80) for reconstrucing the second derivative Radon domain data (R"(?,p)) 
into a three-dimensional Image representation (f(x,y,z)); and 
10 an image memory means (82) for storing the three-dimensional Image representation (f(x,y,2)). 

2. The apparatus as set forth in claim 1 wherein the rebinning means includes a means for mapping from 
preselected points in a cone beam domain to corresponding points in the first derivative Radon domain. 

75 3. The apparatus as set forth In claim 1 wherein the rebinning means includes a means for mapping from 
preselected points in the first derivative Radon domain to corresponding points in a cone beam domain. 

4. The apparatus as set forth in claim 1 including a monitor means for converting a portion of the image 
representation into a human readable display. 

20 

5. The apparatus as set forth in claim 1 further Including a scaling means for scaling the helical data, the 
scaling means being disposed between radiation detection means and the partial derivative means. 

6. The apparatus as set forth In claim 1 wherein the partial derivative means include: 

25 a first partial derivative means for taking a partial derivative along a horizontal direction; 

a second partial derivative means for taking a partial derivative along a vertical direction; and 
a linear combining means for linearly combining the partial derivatives along the horizontal and 
vertical directions to produce the partial derivative along the arbitrary direction. 

30 7. The apparatus as set forth in claim 6 wherein the linear combining means includes a means for 
multiplying the partial derivatives along the vertical and horizontal directions with sine and cosine values 
of an angle between one of the horizontal and vertical directions and the ari^itrary direction in the 
detector plane. 

35 8. The apparatus as set forth in claim 7 wherein for each angle, the integrating means sums row by row to 
produce a one-dimensional array. 

9. A method for generating an image representation of an interior portion of a subject, the method 
comprising: 

40 generating cone beam electrical data from detected radiation travelling along a cone of rays which 

converge at a focal point, the radiation being detected by a radiation detection means along a detection 

plane; 

moving the radiation detecting means such that the focal point moves in a helical orbit; 
taking a partial derivative of the electrical data with respect to an arbitrary direction along the 
45 detection plane; 

projecting the partial derivative along a direction perpendicular to the arbitrary direction along the 
detection plane; 

rebinning the projections Into first derivative Radon domain data; 
taking a second derivative of the Radon domain data; 
50 backprojecting the second derivative Radon domain data into a three-dimensional image repre- 

sentation; 

converting at least a portion of the image representation into a human readable display. 

10. The method as set forth in claim 9 further including storing the three-dimensional image representation. 

55 

11. The method as set forth in claim 9 further including scaling the electrical data prior to taking the partial 
derivative along the arbitrary direction in the detection plane. 
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1i The method as set forth in claim 9 wherein the rebinning includes mapping from preselected points of ' 
the projected partial derivative electrical data along the detection plane to corresponding points in the 
first derivative Radon domain. 

5 13. The method as set forth in claim 9 wherein the rebinning includes mapping from preselected points in 
the first derivative Radon domain to conresponding points on the detection plane. 

14. The method as set forth In claim 9 wherein taking the partial derivative includes: 

taking a partial derivative along a horizontal direction; 
^0 taking a partial derivative along a vertical direction; 

linearly combining the partial derivatives along the horizontal and vertical directions to produce the 
partial derivative along the arbitrary direction. 

15. The method as set forth in claim 14 wherein the linear combining includes multipling the partial 
derivatives along the vertical and horizontal directions with sine and cosine values of an angle between 
one of the horizontal and vertical directions and the arbitrary direction in the detection plane. 

16. The method as set forth in claim 15 wherein for each angle, the projecting step Includes summing row 
by row to produce one-dimensional arrays which are reblnned. 
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